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ABSTRACT 

Wc present new mid-infrared (MIR) imaging data for three Type-1 Scyfcrt galaxies obtained with 
T-ReCS on the Gemini-South Telescope at subarcsecond resolution. Our aim is to enlarge the sample 
studied in a previous work to compare the properties of Type-1 and Type-2 Seyfert tori using clumpy 
torus models and a Bayesian approach to fit the infrared nuclear spectral energy distributions (SEDs). 
Thus, the sample considered here comprises 7 Type-1, 11 Type-2, and 3 intermediate-type Seyferts. 
The unresolved IR emission of the Scyfcrt 1 galaxies can be reproduced by a combination of dust 
heated by the central engine and direct AGN emission, while for the Seyfert 2 nuclei only dust 
emission is considered. These dusty tori have physical sizes smaller than 6 pc radius, as derived from 
our fits. Unification schemes of AGN account for a variety of observational differences in terms of 
viewing geometry. However, we find evidence that strong unification may not hold, and that the 
immediate dusty surroundings of Type-1 and Type-2 Seyfert nuclei are intrinsically different. The 
Type-2 tori studied here are broader, have more clumps, and these clumps have lower optical depths 
than those of Type-1 tori. The larger the covering factor of the torus, the smaller the probability 
of having direct view of the AGN, and vice-versa. In our sample, Seyfert 2 tori have larger covering 
factors (Ct=0.95±0.02) and smaller escape probabilities (Pesc=0.05±g Q3 %) than those of Scyfcrt 1 
(Ct=0.5±0.1; Pesc=18±3 %). All the previous differences are significant according to the KuUback- 
Leibler divergence. Thus, on the basis of the results presented here, the classification of a Seyfert 
galaxy as a Type-1 or Type-2 depends more on the intrinsic properties of the torus rather than on its 
mere inclination towards us, in contradiction with the simplest unification model. 
Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - infrared: galaxies 



1. INTRODUCTION 

Observational evidence in the X-rays and the MIR in- 
dicates that the strong AGN continuum source must be 
absorbed by obscuring material over a wide solid angle 
(see e.g., Antonucci & Miller 1985; Maiolino et al. 1998; 
Risaliti et al. 2002). According to observed spectra of dif- 
ferent AGN types, the obscuring structure has to block 
the emission of the subparsec-scale Broad-Line Region 
(BLR) where the broad lines are produced, but not that 
of the kiloparsec-scale Narrow-Line Region (NLR). 

The unified model for active galaxies (Antonucci 1993; 
Urry & Padovani 1995) is based on the existence of a 
dusty toroidal structure surrounding the central region 
of AGN. This toroidal geometry explains the biconi- 
cal shapes observed in Hubble Space Telescope (HST) 
imaging of several AGNs (Tadhunter & Tsvetanov 1989; 
Malkan et al. 1998; Tadhunter et al. 1999) and also 
the polarimetric observations (Antonucci & Miller 1985; 
Packham et al. 1997). Thus, considering this geometry 
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of the obscuring material, the central engines of Type-1 
AGN can be seen directly, resulting in typical spectra 
with both narrow and broad emission lines, whereas in 
Type-2 AGN the BLR is obscured. 

Pioneering work in modelling the dusty torus (Pier & 
Krolik 1992, 1993; Granato & Danese 1994; Efstathiou & 
Rowan-Robinson 1995; Granato et al. 1997; Siebenmor- 
gen et al. 2004) assumed a uniform dust density distribu- 
tion to simplify the modelling, although from the start, 
Krolik & Begelman (1988) realized that smooth dust dis- 
tributions cannot survive within the AGN vicinity. They 
proposed instead that the material in the torus must be 
distributed in a clumpy structure, in order to prevent the 
dust grains from being destroyed by the hot surrounding 
gas. 

The IR range (and particularly the MIR) is key to set 
constraints on the torus models, since the reprocessed 
radiation from the dust in the torus is re-emitted in this 
range. However, in comparing the predictions of any 
torus model with observations, its small-scale emission 
must be isolated. High angular resolution is then es- 
sential to separate torus emission from stellar emission 
and star-heated dust in the near-IR (NIR) and MIR, re- 
spectively. Indeed, starlight dominates the nuclear NIR 
emission of Seyfert 2 galaxies when using large aperture 
data (see e.g., Alonso-Herrero et al. 1996) and still has 
a significant contribution for Seyfert 1 galaxies (Koti- 
lainen et al. 1992). Similar contamination problems can 
be present in the MIR with the star-heated dust and 
dust in the ionization cones (Alonso-Herrero et al. 2006; 
Mason et al. 2006). 
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Another controversial issue about the torus structure is 
its typical dimensions. Pier & Krolik (1993) and Granato 
& Danese (1994) reproduced the infrared observations of 
nearby Seyfcrt galaxies with ^100 pc scale tori. How- 
ever, hard X-ray observations showed that about half of 
nearby Type-2 Seyferts are Compton-thick (i.e., they are 
obscured by a column density higher than 10^'* cm~^; 
Risaliti et al. 1999). For these highly obscured sources 
the torus dimensions are expected to be of a few parsecs, 
because otherwise the dynamical mass of the obscuring 
material would be too large to be realistic (Risaliti et 
al. 1999). In addition, recent ground-based MIR obser- 
vations of nearby Seyferts reveal that the torus size is 
likely restricted to a few parsecs. Packham et al. (2005) 
and Radomski et al. (2008) established upper limits of 2 
and 1.6 pc for the outer radii of the Circinus galaxy and 
Centaurus A tori, respectively. Besides, interferomet- 
ric observations obtained with the MIR Interferometric 
Instrument (MIDI) at the Very Large Telescope Interfer- 
ometer (VLTI) of Circinus, NGC 1068, and Centaurus A 
suggest a scenario where the torus emission would only 
extend out to R = 1 pc (Tristram et al. 2007), R = 1.7 - 
2 pc (Jaffe et al. 2004; Raban et al. 2009), and R = 0.3 
pc (Meisenheimer et al. 2007), respectively. 

In order to solve the discrepancies between observa- 
tions and previous models, an intensive search for an al- 
ternative torus geometry has been carried out in the last 
decade. The first results of radiative transfer calculations 
of a clumpy rather than a smooth medium were reported 
by Nenkova et al. (2002) and Elitzur & Shlosman (2006), 
and further work was done by DuUcmond & van Bemmel 
(2005). The clumpy dusty torus models (Nenkova et al. 
2002, 2008a,b; Honig et al. 2006; Schartmann et al. 2008) 
propose that the dust is distributed in clumps, instead 
of homogeneously filling the torus volume. These mod- 
els are making significant progress in accounting for the 
MIR emission of AGNs (Mason et al. 2006, 2009; Mor 
et al. 2009; Horst et al. 2008, 2009; Nikutta et al. 2009; 
Ramos Almeida et al. 2009a; Honig et al. 2010). 

In our previous work (Ramos Almeida et al. 2009a; 
hereafter RA09), we constructed subarcsecond resolu- 
tion IR SEDs for eighteen Seyfert galaxies, mostly Type- 
2 Seyferts. From the comparison between our high an- 
gular resolution MIR fluxes and large aperture data, 
such as those from ISO, IRAS, or Spitzer, we confirmed 
that the former provide a spectral shape that is sub- 
stantially different from that of the large aperture data 
(Rodriguez Espinosa & Perez Garcia 1997). Since our 
nuclear measurements allowed us to better characterize 
the torus emission, we modelled our SEDs with clumpy 
torus models. In general, we found that Type-2 views 
are more inclined than those of Type-Is, and more im- 
portantly, we derive larger covering factors for the Type- 
2 tori (i.e., more clumps and wider torus angular dis- 
tributions). This would imply that the observed differ- 
ences between Type-1 and Type-2 AGN would not be due 
to orientation effects only, but to intrinsic differences in 
their tori. However, due to the limited size of the sample 
analyzed by RA09, and in particular of Type-1 Seyferts 
compared with Type-2s, our aim is to enlarge the sample 
studied in the previous work with new Seyfert 1 infrared 
data to compare the properties of Type-1 and Type-2 
Seyfert tori. 

In this work, we report new subarcsecond MIR imaging 



data for the 3 nearby Type-1 Seyfert galaxies NGC 7469, 
NGC 6221, and NGC 6814, for which we estimate unre- 
solved nuclear MIR fluxes. We enlarge the sample by in- 
cluding the galaxies NGC 1097, NGC 1566, NGC 3227, 
and NGC 4151, which have similar MIR data, and we 
compile NIR nuclear fluxes from the literature of similar 
resolution to construct nuclear SEDs for all the galaxies. 
We fit these SEDs with clumpy torus models which we 
interpolate from the Nenkova et al. (2008a, b) database, 
and compare them with the larger sample studied by 
RA09. Table 1 summarizes key observational properties 
of the sources in the sample. Section 2 describes the ob- 
servations, data reduction, and compilation of NIR and 
MIR fluxes. Sections 3 and 4 present the main obser- 
vational results, and in §5 we report the modelling re- 
sults. We discuss differences between Type-1 and Type- 
2 Seyferts and draw conclusions about the clumpy torus 
models and AGN obscuration in general in §6. Finally, 
Section 7 summarizes the main conclusions of this work. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. MIR Imaging Observations 

In order to enlarge the sample of 18 Seyfert galax- 
ies presented in RA09, which comprises 12 Seyfert 2 
(Sy2), two Seyfert 1.9 (Syl.9), one Seyfert 1.8 (Syl.8), 
two Seyfert 1.5 (Syl.5), and one Seyfert 1 galaxy (Syl), 
we obtained new subarcsecond MIR observations of the 
Type-1 Seyferts NGC 6221, NGC 6814, and NGC 7469 
(see Table 1). 

The observations were performed with the MIR 
camera/spectrograph T-ReCS (Thermal- Region Camera 
Spectrograph; Telesco et al. 1998) on the Gemini-South 
telescope during the Summer of 2009. T-ReCS uses a 
Raytheon 320x240 pixel Si:As IBC array, providing a 
plate scale of 0.089" pixel" ^, corresponding to a FOV 
of 28.5"x21.4". The filters employed for the observations 
were the narrow Si-2 filter (Ac=8.74 //m, AA=0.78 ^m, 
50% cut-on/off) and the broad Qa filter (Ac=18.3 ^m, 
AA=1.5 /xm, 50% cut-on/off). The resolutions obtained 
were 0.3" at 8.7 /im and 0.5" at 18.3 /im, as measured 
from the observed Point Spread Function (PSF) stars. A 
summary of the observations is reported in Table 2. 

The standard chopping-nodding technique was used to 
remove the time-variable sky background, the telescope 
thermal emission, and the so-called 1/f detector noise. 
The chopping throw was 15", and the telescope was nod- 
ded 15" in the direction of the chop every 45 s. The 
difference for each chopped pair for each given nodding 
set was calculated, and the nod sets were then differenced 
and combined until a single image was created. The data 
were reduced using in-house-developed IDL routines. 

Observations of flux standard stars were made for the 
flux calibration of each galaxy through the same filters. 
The uncertainties in the flux calibration are typically '^5- 
10% at 8.7 /im and ~15-20% at 18.3 fim. PSF star 
observations were also made immediately prior to or af- 
ter each galaxy observation to accurately sample the im- 
age quality. These images were employed to determine 
the unresolved (i.e., nuclear) component of each galaxy. 
The PSF star, scaled to the peak of the galaxy emis- 
sion, represents the maximum contribution of the unre- 
solved source (100%), where we integrate flux within an 
aperture of 2". The residual of the total emission minus 
the scaled PSF represents the host galaxy contribution, 
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TABLE 1 
Basic Galaxy Data 



Galaxy 


Seyfert Type 


Ref. 


z 


Distance 


Scale 


Ref. 










(Mpc) 


(pc arcscc^^) 




NGC 1097 


Sylt 


Al 


0.0042 


19 


92 


Bl 


NGC 1566 


Syl 


A2 


0.0050 


20 


97 


B2 


NGC 6221 


Syl 


A3 


0.0050 


18 


87 


B3 


NGC 6814 


HII/Syl.5 


A4 


0.0052 


21 


102 


B4 


NGC 7469 


Syl 


A5 


0.0163 


65 


315 


B5 


NGC 3227 


Syl. 5 


A6 


0.0039 


17 


82 


B6 


NGC 4151 


Syl.5 


A7 


0.0033 


13 


64 


B7 



References. — (Al) Storchi-Borgmann et al. (1997); (A2) Kriss et al. 
(1990); (A3) Levenson et al. (2001); (A4) Veron-Cetty & Veron (2006); (A5) 
Osterbrock & Martel (1993); (A6) Rubin & Ford (1968); (A7) Ayani & Mae- 
hara (1991); (Bl) Willick et al. (1997); (B2) Sandage & Bedke (1994); (B3) 
Koribalski & Dickey (2004); (B4) Liszt & Dickey (1995); (B5) Heckman et al. 
(1986); (B6) Garcia (1993); (B7) Radomski et al. (2003). 

Note. — Classification and distance arc taken from the literature (refer- 
ences below) and spectroscopic redshift from the NASA/IPAC Extragalactic 
Database (NED). 

tOriginally classified as a LINER by Keel (1983) and Phillips et al. (1984). 



TABLE 2 
Summary of MIR Observations 



Galaxy Filters Observation On-Source Time (s) PSF FWHM 

epoch N band Q band N band Q band 



NGC 1097 


Si-5. ( 


3a 


Sep 


2005 


456 


912 


0.41" 


0.52" 


NGC 1566 


Si-2, ( 


3a 


Sep 


2005 


152 


304 


0.30" 


0.53" 


NGC 6221 


Si-2, ( 


3a 


Aug 


2009 


145 


203 


0.32" 


0.55" 


NGC 6814 


Si-2, ( 


3a 


Aug 


2009 


145 


203 


0.28" 


0.53" 


NGC 7469 


Si-2, ( 


3a 


Sep 


2009 


145 


203 


0.31" 


0.55" 


NGC 3227 


N' 




Apr 


2006 


300 




0.39" 




NGC 4151 


N, IHW18 


May 


2001 


360 


480 


0.53" 


0.58" 



Note. — Images were obtained in the 8.74 /^m (Si-2, AA = 0.78 /im at 50% 
cut-on/off), 11.66 ^m (Si-5, AA = 1.13 /im), and 18.30 ^m (Qa, AA = 1.5 /^m) 
T-ReCS filters. NGC 3227 was observed in the 11.29 (N', AA = 2.4 /^m) 
Michelle/Gemini-North filter, and NGC 4151 in the 10.75 /^m (N, AA =5.2 ^im) 
and 18.17 /xm (IHW18, AA =1.7 ^tm) OSCIR/Gemini-North filters. 



which is analyzed in Section 3. We require a flat profile in 
the residual for a realistic galaxy profile over the central 
pixels. Therefore, we reduce the scale of the PSF from 
matching the peak of the galaxy emission (100%), when 
necessary, to obtain the unresolved fluxes reported in Ta- 
ble 3. They include corresponding aperture corrections 
to take into account possible flux losses when integrating 
the scaled PSF flux in the 2" aperture. 

Figure 1 shows an example of PSF subtraction at var- 
ious levels (in contours) for NGC 7469 in the Si-2 T- 
ReCS filter, following Radomski et al. (2002, 2003). The 
residual profiles from the different scales demonstrate the 
best-fitting result. The uncertainty in the unresolved 
fluxes determination from PSF subtraction is '^10-15%. 
Thus, we estimated the errors in the flux densities re- 
ported in Table 3 by adding quadratically the flux cal- 
ibration and PSF subtraction uncertainties, resulting in 
-15% at 8.7 ^m and ~25% at 18.3 \mi. 

2.2. Compilation of NIR and MIR High Spatial 
Resolution Data 

In addition to the new T-ReCS observations described 
in Section 2.1, we include here two more Syl galaxies 
with aheady pubhshcd T-ReCS data: NGC 1097 and 



TABLE 3 
Unresolved MIR Fluxes 



Galaxy Level of PSF subtraction Flux Density (mjy) 





N band 


Qa band 


N band 


Qa band 


NGC 1097 


80% 


20% 


23 


55 


NGC 1566=^ 


100% 


100% 


29 


117 


NGC 6221 


80% 


30% 


48 


314 


NGC 6814 


100% 


40% 


53 


159 


NGC 7469 


85% 


100% 


174 


1354 


NGC 3227 


100% 




401 




NGC 4151 


90% 


100% 


1320 


3200 



Note. — The percentages of PSF subtraction level are re- 
ported in the employed filters (listed in Table 2). Errors in flux 
densities are dominated in general by uncertainties in the flux 
calibration and PSF subtraction (~15% at N band and ~25% 
at Qa band). Fluxes include aperture corrections. 

''MIR fluxes for NGC 1566 are from aperture photometry using 
0.45" aperture radii ( RA09). 

NGC 1566^. We reduced the images presented in Ma- 
son et al. (2007) for NGC 1097 and obtained unresolved 

3 NGC 1097 and NGC 1566 were not included in the analysis of 
the SEDs presented in RA09 because of the lack of high angular 
resolution NIR fluxes for them at the time of publication. 
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Fig. 1. — 8.74 /im contour plots of NGC 7469 at the 3cr level, the PSF star, and scaled subtraction of the PSF for this galaxy at various 
levels (100%, 90%, and 80%). The residuals of the subtraction in the lower right panel show the host galaxy profile. 



MIR fluxes by performing the same technique explained 
in Section 2.1 (see Table 3). For NGC 1566 we compiled 
the MIR nuclear fluxes reported RA09. Both galaxies 
were observed in September 2005 with T-ReCS: NGC 
1566 in the Si-2 and Qa filters, and NGC 1097 in the 
Si-5 (Ac=11.66 /im, AA=1.13 /tm, 50% cut-on/off) and 
Qa filters (see Table 2). The resolutions of the images 
are 0.3" at 8.7 ^lm, ~0.4" at 11.66 /im, and 0.5" at 18.3 
/zm. 

NIR subarcsecond resolution nuclear fluxes compiled 
from the literature are reported in Table 4. For NGC 
1097, NGC 1566, and NGC 7469, Prieto et al. (2010) 
reported diffraction-limited and near-diffraction-limited 
adaptive optics NACO/VLT fluxes. The three galaxies 
are unresolved in the NIR down to the highest resolution 
achieved (FWHM-0.15" for NGC 1097 in the L'-band, 
-0.12" for NGC 1566 in the K-band, and -0.08" for 
NGC 7469 in the H-band). 

For NGC 7469, Prieto et al. (2010) reported NACO 
J-, H-, and K-band nuclear fluxes obtained from obser- 
vations in Nov 2002, as well as in the L' and NB-4.05 
/im filters, observed in this case in Dec 2005. First, we 
discarded the narrow-band filter, since it is designed to 
collect the Bra emission, which in the case of this galaxy 
is important (as inferred from the Br7 line detected with 
NIR spectroscopy; Genzel et al. 1995; Hicks & Malkan 
2008). The L'-band filter also contains Bra, which is very 
likely compromising the flux measurement, since the L' 
and NB-4.05 fluxes do not match the SED shape of the 
remaining J, H, K, Si-2 and Qa data. Indeed, Prieto et 
al. (2010) also reported a NICMOS/HST flux measure- 
ment in the fiher F187N, obtained in 2007. This data 
point lies in between the NACO H and K measurements 
in wavelength and flux, which gives us extra-confidence 



in the reliability of the NACO J, H, and K measure- 
ments. Considering all the previous, we finally decided 
to consider the NACO L'-band flux as an upper limit^*^. 

For NGC 6221 and NGC 6814, which do not have 
any published subarcsecond resolution NIR data, we 
retrieved broad-band images from the Hubble Legacy 
Archive" obtained with NICMOS on the HST. The 
two galaxies were observed in the F160W filter, using 
the NIC2 camera, as part of the program 7330 (PI: 
Mulchaey, J.). The typical FWHM for an unresolved 
PSF is -0.13" using the F160W filter with NIC2. De- 
tails of the observations can be found in Scoville et al. 
(2000) and Regan & Mulchaey (1999). For the anal- 
ysis, we first separated the nuclear emission from the 
underlying host galaxy emission. We then applied the 
two-dimensional image decomposition GALFIT program 
(Peng et al. 2002) to fit and subtract the unresolved com- 
ponent (PSF). PSF models were created using the Tiny- 
Tim software package, which includes the optics of HST 
plus the specifics of the camera and filter system (Krist 
1993). We checked that no prominent emission lines are 
included in the wavelength range covered by the filter 
F160W. The resuhing unresolved NIR fluxes for NGC 
6221 and NGC 6814 are reported in Table 4. 

We finally include the Syl.5 galaxies NGC 3227 and 
NGC 4151 in this study, which were also part of the 
RA09 sample. The MIR nuclear fluxes are the same re- 
ported in the latter work (see description of the observa- 

^'^ The nucleus of NGC 7469 had undergone different periods 
of activity, including a maximum in the optical happening in the 
period 1996 to 2000, followed by a relaxation epoch in the following 
years (Prieto et al. 2010). For this reason, here we only consider 
data obtained from 2002. 

http;//archivc. stsci.edu/ 
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tions in Section 2.1 in RA09), whereas the NIR fluxes 
have been updated (see Table 4). 

Using the NIR nuclear fluxes reported in Table 4 in 
combination with our MIR unresolved measurements 
(Table 3) we construct AGN-dominated SEDs for the 
iivc Syl galaxies. 

3. THE MIR EXTENDED EMISSION OF NGC 7469 AND 
NGC 6221 

Our new MIR imaging data reveal complex extended 
emission for NGC 7469 and NGC 6221, which is known 
to be intense and associated with emitting-dust heated 
by star formation. On the other hand, NGC 6814 lacks 
of any extended emission. In this section we present the 
MIR images of NGC 7469 and NGC 6221, and compare 
them with published data in different wavelength ranges. 

3.1. NGC 7469 

The most spectacular morphological feature of this Syl 
galaxy is a circumnuclear ring of powerful starbursts of 
~1.6 kpc diameter, which is deeply embedded in a large 
cloud of molecular gas and dust. This ring contains sev- 
eral super star clusters and regions of star formation and 
it accounts for two-thirds of the galaxy bolometric lu- 
minosity (Genzel et al. 1995). The star-forming ring 
has been the subject of study of several works in differ- 
ent wavelengths (see Dfaz-Santos et al. 2007, hereafter 
DS07, and references therein), including the MIR (Miles 
et al. 1994; Soifer et al. 2003; Horst et al. 2009). Based 
on VISIR/VLT MIR observations at 12.3 ^m, Horst et 
al. (2009) report the detection of distinct knots of star 
formation around the nucleus located at a distance of 
~1.3" (^^400 pc), although with low signal-to-noise. Us- 
ing NIR HST data, DS07 identified 30 clusters of star 
formation in the ring at 1.1 fim. By fitting the individual 
ultraviolet-to-NIR SEDs of the clusters, the authors re- 
ported the presence of a dominant intermediate age pop- 
ulation (8-20 Myr) and a younger and more extinguished 
one (■~l-3 Myr). The latter does not coincide with the 
optical/NIR continuum-emitting regions, but seems to 
be traced by the MIR/radio emission, less affected by 
extinction than the optical/NIR. 

Figure 2 shows the high resolution 8.74 and 18.3 ^m T- 
ReCS images of NGC 7469. At both wavelengths we 
detect extended emission with high signal-to-noise coin- 
cident with the star-forming ring. Indeed, our flux maps 
appears very similar to the high resolution 11.7 fj,m con- 
tours presented by Miles et al. (1994) and to the high 
resolution (0.2") VLA 8.4 GHz radio map presented in 
Colina et al. (2001). We identified the brightest knots in 
our MIR images in Figure 2 using the same notation as 
in Miles et al. (1994), where the AGN is labelled A. The 
B and C knots in our images correspond to the bright- 
est regions in radio, according to the 5 GHz (Wilson et 
al. 1991) and 8.4 GHz radio maps (Colina et al. 2001). 
We also identified the knots D, E, and F. In Table 5 we 
report the star clusters identified by DS07 which better 
match the positions of the A to F MIR compact regions. 
The distances between these knots and the AGN range 
from 1.4" to 1.8" (median distance of ~480 pc). Ah 
the knots appear more compact in the 18.3 fim image 
than in the 8.74 //m, where the ring emission is more ex- 
tended. This is expected since the 8.74 fim filter contains 
the 8.6 /im Polycyclic Aromatic Hydrocarbon (PAH) fea- 



ture, whereas the 18.3 fim filter is mostly probing hot 
dust emission. As found by Spitzer, the ~8 fim PAH 
emission of nearby galaxies appears to be more extended 
than at ^^24 fim (e.g., Hclou et al. 2004; Calzetti et al. 
2005). 

In Table 5 we report aperture fluxes for the six iden- 
tified knots in both the 8.74 and 18.3 fim images cal- 
culated using IRAF^^. The aperture radius was defined 
on the basis of the resolution element in the Qa band 
(0.55"), and corresponding aperture corrections were ap- 
plied, since the knots are only barely resolved in both 
bands. Positions relative to the nucleus (A) in arcsec- 
onds are also given. 

3.2. NGC 6221 

The galaxy NGC 6221 is a prototypical example of 
the so-called "X-ray-loud composite galaxies" (see e.g., 
Moran et al. 1996). This classification comes from 
the comparison between its optical spectrum, which is 
starburst-like, and its X-ray data, where the AGN is re- 
vealed (Levenson et al. 2001). According to the X-ray 
data, and in particular the width of the Fe Ka line, the 
orientation of the Seyfert nucleus is Type-1 (Levenson 
et al. 2001). However, as mentioned above, the optical 
spectrum of the galaxy resembles more that of a starburst 
galaxy, implying that there must be a big amount of dust 
(likely associated with the starburst) along the line of 
sight (LOS) hiding the BLR. A nuclear optical extinction 
of Ay =3.0 mag was measured from the optical spectrum 
by Levenson et al. (2001). They also presented HST im- 
ages of the central region of NGC 6221 in the optical 
(F606W/WFPC2) and in the NIR (F160W/NICMOS) 
and identified the bright central NIR source as the AGN. 
At optical wavelengths, the nucleus is diffuse and weaker 
than other bright knots, identified as star clusters. 

Our 8.7 and 18.3 fim images of NGC 6221 are shown 
in Figure 3. They reveal spectacular extended emission 
with two bright knots. The AGN is centered in both 
images. The other bright knot, at ^1.9" SW from the 
nucleus, which is roughly coincident with the SW star- 
burst region shown in the F606W optical image in Figure 

3 of Levenson et al. (2001), reaches practically the same 
intensity as the AGN at 8.7 fim and appears brighter at 
18.3 fim. We measured aperture fluxes in a 0.7" radius 
for the SW knot, selected to collect the bulk of its MIR 
emission, and obtain fluxes of 42 mJy at 8.7 fim and 422 
mjy at 18.3 fim. Surrounding the AGN there is more 
diffuse emission extending towards the North, which in 
this case is more intense in the 8.7 fim image than in the 
18.3 ^m one. As already mentioned in Section 3.1, the 
~8 fim PAH emission of nearby galaxies appears more 
extended than the ~24 /im emission (e.g., Helou et al. 
2004; Calzetti et al. 2005). 

4. SED OBSERVATIONAL PROPERTIES 

4.1. Average Seyfert 1 Spectral Energy Distribution 

Using the MIR and NIR data reported in Tables 3 and 

4 we construct subarcsecond resolution nuclear SEDs in 
the wavelength range from ^1 to 18 fim for the five Type- 
1 Seyferts analyzed here. Figure 4 shows a comparison 

^■^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
the Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (http://iraf.noao.edu/). 
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TABLE 4 

High Spatial Resolution NIR Fluxes 



Galaxy 


J band 


Flux Density (mjy) 
H band K band 


L band 


Filters 


Reference(s) 


NGC 1097 


l.litO.l 


2.7±0.1 


3.9±0.1 


llitl 


NACO J,H,K,L' 


a 


NGC 1566 


l.litO.l 




2.1±0.1 


7.8±0.1 


NACO J,K,L' 


a 


NGC 6221 




2.1±0.2 






F160W 


b 


NGC 6814 




5.2±0.5 






F160W 


b 


NGC 7469=^ 


S.OitO.l 


15±1 


20±1 


84±1 


NACO J,H,K,L' 


a 


NGC 32271^ 
NGC 4151'^ 


60±6 


7.8±0.8 
lOOitlO 


16.4±1.7 
197±20 


46.7±9.3 
325±65 


F160W,F222M, NSFCam L 
F110W,F160W,F222M, NSFCam L 


c,d 

c,e 



References. — (a) Prieto et al. (2010); (b) This work; (c) Kishimoto ct al. (2007); (d) Alonso-Herrero et 
al. (2003); (e) Ward et al. (1987) 

Note. — Ground-based instruments and telescopes are NACO on the 8 m VLT and NSFCam on the 3 m 
NASA IRTF. Measurements in the FllOW, F160W, and F222M filters are from NICMOS on HST. 

''Prieto et al. (2010) also reported a flux measurement of 19±1 mJy obtained from a NICMOS/HST observation 
in the filter F187N. 

''Ward et al. (1987) also reported M-band flux measurements of 72±27 mJy for NGC 3227 and 449±34 mJy 
for NGC 4151 obtained with IRCAM3 on the 3.8 m UKIRT. 




Fig. 2. — 8.74 fim (left) and 18.3 /im (right) T-ReCS images of NGC 7469 smoothed by using a moving box of 3 pixel size. The images 
size is 7.8 arcsec side. North is up, and East to the left. Emitting regions identified in both the Si-2 and Qa images, and also coincident 
with the radio emission, arc labelled from A to F, where A corresponds to the AGN. 

between their spectral shapes and the average Sy2 SED 
from RA09. This mean template was constructed using 
individual Sy2 data of the same angular resolution as 
that achieved in this work (<0.55"). In the same way, 
we have constructed an average Type-1 Seyfert template 
using the IR nuclear SEDs of the seven galaxies studied 
here^^. The spectral shape of Syl and Syl.5 galaxies 
is practically identical (as can be seen from Figure 4). 
Based on the previous, and on the fact that both types 
of nuclei present broad lines in their optical spectra, in 
this work wc consider them as Typc-1 Seyferts. 

The Sy2 template defines the wavelength grid, and we 
performed a quadratic interpolation of nearby measure- 
ments of the individual Syl galaxies onto its scale (1.265, 

^'^ We have not considered the NIR ground-based data reported 
in Table 4 of RA09 for NGC 3227 and NGC 4151 in the con- 
struction of the mean template to be consistent with the angular 
resolutions of the other SEDs. 



1.60, 2.18, 3.80, 4.80, 8.74, and 18.3 ^m). We did not in- 
terpolate the sparse observations of NGC 6221 and NGC 
6814, which only have NICMOS 1.6 fim data in addi- 
tion to the MIR measurements. The interpolated fluxes 
were used solely for the purpose of deriving the average 
Syl template. The error bars correspond to the stan- 
dard deviation of each averaged point, except for the 
8.74 fim point (the wavelength chosen for the normal- 
ization). In this case, we assigned a 15% error, which is 
the nominal percentage considered for the N-band flux 
measurements (see Section 2). 

We measured the 1.265-18.3 /im IR slope (fy a v'"^'" ) 
of the Syl template, which is representative of the indi- 
vidual Syl SEDs: am = 1.7±0.3. We also calculated the 
NIR {aNiB. = 1.6±0.2, from 1.265 to 8.74 fim) and MIR 
spectral indexes {aMiR = 2.0 ± 0.2. using the 8.7 and 
18.3 yum points). A flat NIR slope indicates an impor- 
tant contribution of hot dust emission (up to ^1000-1200 
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TABLE 5 

Positions and flux densities of NGC 7469 MIR knots 



Knot 


ID 




8.74 


fim 




18.3 






DS07 


(X) 


(Y) 


Flux (mjy) 


(X) 


(Y) 


Flux (mJy) 


A 




0.00 


0.00 


229 


0.00 


0.00 


1320 


B 


C20 


-1.58" 


-0.89" 


11.9 


-1.54" 


-0.93" 


50 


C 


C6 


1.09" 


1.07" 


13.4 


0.95" 


1.06" 


79 


D 


CIO 


-0.49" 


1.31" 


11.5 


-0.58" 


1.38" 


107 


E 


C7,C15 


-0.40" 


-1.36" 


15.8 


-0.52" 


-1.49" 


104 


F 


C12 


1.25" 


-0.89" 


3.7 


1.47" 


-0.83" 


14.8 



Note. — Columns (1) and (2) indicate the label asigned to each knot in Figure 
2 and in DS07. B, C, and D are also coincident with the R3, R2, and Rl regions 
in the VLA 8.4 GHz radio maps in Colina ct al. (2001). Columns (3,4) and (6,7) 
indicate the position of each knot relative to A at 8.7 and 18.3 ^m, respectively. 
Columns (5) and (8) list the flux densities of the knots in a 0.55" aperture radius. 
Fluxes include aperture corrections of 19% at 8.7 fim and 37% at 18.3 fim. Errors 
in flux densities are dominated by uncertainties in the flux calibration (~5-10% at 
8.7 ^tm and ~15-20% at 18.3 fim). 




Fig. 3. — 8.74 fim (left) and 18.3 ^m (right) T-ReCS images of NGC 6221 smoothed by using a moving box of 3 pixel size. The images 
size is 7.8 arcsec side. North is up, and Elast to the left. The AGN is centered in both images. 



K; Rieke & Lebofsky 1981; Barvainis 1987) that comes 
from the immediate vicinity of the AGN. am, ctNiR, and 
otMiR values for the individual Type-1 Seyfert galaxies 
and for the mean Syl and Sy2 SEDs are shown in Ta- 
ble 6. The shape of the Sy2 mean SED is very steep 
{aiB, = 3.1±0.9, aNiR = 3.6±0.8, and umir = 2.0±0.2) 
compared with those of the Type-1 Seyferts. In general 
Sy2 have steeper 1-10 ^m SEDs than Syl (Rieke 1978 
Edelson et al. 1987; Ward et al. 1987; Fadda et al. 1998 
Alonso-Herrero et al. 2001, 2003). On the contrary, the 
MIR slope results to be the same (aA//fl = 2.0 ± 0.2) for 
both the Syl and Sy2 templates. 

Alonso-Herrero et al. (2003) also reported IR spectral 
indices measured from 1 to 16 /Ltm for Syl and Syl. 5 
galaxies (a^j^°"'*°=l. 5-1.6). On the other hand, the NIR 
slopes of the Syl. 8 and Syl. 9 in RA09 have interme- 
diate values between those of Sy2 and Syl (mean slope 
aiH = 2.0 ± 0.4), also in agreement with the IR slopes 
reported in Alonso-Herrero et al. (2003) for Syl. 8 and 
Syl.9 (a^Jj°""°=1.8-2.6). In summary, the slope of the 



IR nuclear SED is generally correlated with Seyfert type: 
Type-2 nuclei show steeper SEDs, whereas Type-1 and 
intermediate Seyferts are flatter. The NIR excess respon- 
sible for flattening the SED of the Type-1 nuclei would 
come from the contribution of hot dust in the directly- 
illuminated faces of the clouds in the torus, as well as 
from the direct AGN emission (i.e., the tail of the opti- 
cal power-law continuum). 

A similar comparison between Type-1 and Type-2 
spectral shapes can be done by using the H/N and N/Q 
flux ratios. H/N is larger for Type-1 Seyferts (0.07±0.03) 
than for Sy2 (0.003±0.002), as measured from the indi- 
vidual values of the Syl considered here and the Sy2 
in RA09. The difference in H/N between Syl and Sy2 
galaxies is significantly different at the 100% confidence 
level, according to the Kolmogorov-Smirnov (K-S) test. 
This ratio depends on both the torus inclination and cov- 
ering factor, as we will discuss in Section 6.1. On the 
other hand, N/Q is very similar for Type-1 and Type-2 
Seyferts (mean values of 0.27±0.11 and 0.23±0.14 respec- 
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TABLE 6 
Spectral Shape Information 



Galaxy 




OlNIR 


OlMIR 


H/N 


N/Q 


Average Sy2 


3.1±0.9 


3.6±0.8 


2.0±0.2 


0.003 


0.23 


Average Syl 


1.7±0.3 


1.6±0.2 


2.0±0.2 


0.071 


0.22 


NGC 1097 


1.4±0.2 


1.6±0.2 


1.2±0.1 


0.117 


0.42 


NGC 1566 


1.8±0.3 


1.8±0.3 


1.9±0.2 


0.040=' 


0.25 


NGC 3227 


2.0±0.4 


2.0±0.4 




0.019 




NGC 4151 


1.4±0.2 


1.4±0.1 


1.5±0.1 


0.076 


0.41 


NGC 6221 


2.0±0.5 


1.8±0.4 


2.5±0.3 


0.044 


0.15 


NGC 6814 


1.4±0.2 


1.4±0.2 


1.5±0.1 


0.098 


0.33 


NGC 7469 


1.8±0.3 


1.5±0.2 


2.8±0.4 


0.086 


0.13 



Fig. 4. — Observed IR SEDs for the seven Typc-1 Scyfert galaxies (in color and with different symbols) used for the construction of 
the Syl average template (dashed pink). The Sy2 template from RA09 (solid black) is also plotted for comparison. The SEDs have been 
normalized at 8.74 fim, and the average Syl and Sy2 SEDs have been shifted in the Y-axis for clarity. 

tributed with a radial extent Y = Ro/ Rd, where Ro and 
Rd are the outer and inner radius of the toroidal distribu- 
tion, respectively (see Figure 5). The inner radius is de- 
fined by the dust sublimation temperature {Tsut ~ 1500 
K), with Rd = 0.4 (1500 K T-l)'^-^{L/10'^^ erg s-'^f-^ 
pc. Within this geometry, each clump has the same op- 
tical depth (ry, defined at the y-band). The average 
number of clouds along a radial equatorial ray is Nq. 
The radial density profile is a power-law (oc r^'^). A 
width parameter, cr, characterizes the angular distribu- 
tion of the clouds, which has a smooth edge. The num- 
ber of clouds along the LOS at an inclination angle i is 
A^LOs(j) = No e(-('-3°)'/'^'). Finally, the optical extinc- 
tion produced by the torus along the LOS is computed 
as A^O^ = 1.086 No ry e(-(*-90)V'^') mag. For a de- 
tailed description of the clumpy models see Nenkova et 
al. (2002, 2008a,b). 

The clumpy database now contains 1.2x10^ models, 
calculated for a fine grid of model parameters. The in- 
herent degeneracy between these parameters has to be 
taken into account when fitting the observables. To this 
end, we recently developed a Bayesian inference tool 
(BayesClumpy) , that extracts as much information as 
possible from the observations. Details on the interpo- 
lation methods and algorithms employed can be found 
in Asensio Ramos & Ramos Almeida (2009). Thus, for 
the following analysis of the Seyfert SEDs, we are not us- 
ing the original set of models described in Nenkova et al. 
(2008a, b), but an interpolated version of them (See Fig- 
ures 3 and 4 in Asensio Ramos & Ramos Almeida 2009 
for a comparison between the original and interpolated 



Note. — Columns 2, 3, and 4 show the fitted spectral 
indexes (f^ a u~") measured from the average Sy2 and Syl 
templates and for the individual Syl SEDs in the whole range 
{am, from ~1 to 18 /^m), in the NIR (ajv/ij, from ~1 to ~9 
lim), and in the MIR {cej^uji, using the N and Q band data 
points) respectively. Columns 5 and 6 give the values of the 
H/N and N/Q band ratios. 

"Due to the lack of H-band data for NGC 1566, we used the 
interpolated value at 1.6 /im for calculating the H/N. 

tively). Values of H/N and N/Q for the individual Syl 
galaxies and the average templates are reported in Table 
6. 



5. SED MODELLING 



5.1. 



Clumpy Dusty Torus Models and Bayesian 
approach 

The clumpy dusty torus models of Nenkova et al. 
(2002) hold that the dust surrounding the central engine 
of an AGN is distributed in clumps, instead of homoge- 
neously filling the torus volume. These clumps are dis- 
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Fig. 5. — Scheme of the clumpy torus described in Nenkova et 
al. (2008a, b). The radial extent of the torus is defined by the outer 
radius (Ro) and the dust sublimation radius (Rd)- All the clouds 
are supposed to have the same Ty, and cr characterizes the width 
of the angular distribution. The number of cloud encounters is 
function of the viewing angle, i. 

models). In this work we use the most up-to-date ver- 
sion of the models that are corrected for a mistake in the 
torus emission calculations, which in principle only af- 
fects the AGN scaling factor (see erratum by Nenkova et 
al. 2010). The present version of BayesClumpy has also 
been updated to use the Multinest algorithm of Feroz 
et al. (2009) for sampling. This algorithm is very ro- 
bust and efficient when sampling from complex posterior 
distributions. 

The prior distributions for the model parameters are 
assumed to be truncated uniform distributions in the in- 
tervals reported in Table 7. Therefore, we give the same 
weight to all the values in each interval. Apart from the 
six parameters that characterize the models, there is an 
additional parameter that accounts for the vertical dis- 
placement required to match the fluxes of a chosen model 
to an observed SED, which we allow to vary freely. This 
vertical shift scales with the AGN bolometric luminosity 
(see Section 6). In order to compare with the observa- 
tions, BayesClumpy simulates the effect of the employed 
filters on the simulated SED by integrating the product 
of the synthetic SED and the filter transmission curve. 
Observational errors are assumed to be Gaussian or up- 
per/lower limit detections. A detailed description of the 
Bayesian inference applied to the clumpy models can be 
found in Asensio Ramos & Ramos Almeida (2009). Ad- 
ditionally, to see an example of the use of clumpy model 
fitting to IR SEDs using BayesClumpy, see RA09. 



5.2. Model Results 



5.2.1. 



Seyfert 1 Individual Fits 

The results of the fitting process of the IR SEDs with 
the interpolated version of the clumpy models of Nenkova 
et al. (2008a, b) are the posterior distributions for the six 
free parameters that describe the models and the verti- 
cal shift. These are indeed the probability distributions 
of each parameter, represented as histograms. When the 
observed data introduce sufficient information into the 
fit, the resulting posteriors will clearly differ from the 
input uniform priors, either showing trends or being cen- 



tered at certain values within the intervals considered. 
For all the Syl fits, we considered uniform priors in the 
intervals shown in Table 7. The only exception is NGC 
7469, for which we use a gaussian prior of 85°±2°for the 
inclination angle of the torus, based on the value of the 
accretion disk viewing angle deduced from X-ray obser- 
vations (Nandra et al. 2007), assuming that the disk and 
the torus are coplanar. 

We fit the individual Syl SEDs with BayesClumpy, 
modelling the torus emission and the direct AGN contri- 
bution (the latter as a broken power law). We also con- 
sider the IR extinction curve of Chiar & Tielens (2006) 
to take into account any possible foreground extinction 
from the host galaxy. The AGN scales self-consistently 
with the torus flux, and the foreground extinction (sep- 
arate from the clumpy torus) is another free parameter, 
which we set as a uniform prior ranging from Ay^O to 
10 mag^^. 

In addition to the Gemini MIR unresolved fluxes re- 
ported in Table 3 we consider MIR nuclear fluxes from 
VISIR compiled from the literature when available (see 
Table Al in Appendix A) We flnd good agreement be- 
tween the VISIR and T-ReCS unresolved fluxes. We did 
not include measurements in the VISIR PAH fllters (8.59, 
11.25 and 11.88 ^i^Jl) for the galaxies NGC 7469 and NGC 
1097 because of their intense star formation and their al- 
ready well-sampled SED. 

Although the solutions to the Bayesian inference prob- 
lem are the probability distributions of each parameter, 
we can translate the results into corresponding spectra 
(Figures 6 and 7). The solid lines correspond to the 
model described by the combination of parameters that 
maximizes their probability distributions (maximum-a- 
posteriori; MAP). Dashed lines represent the model com- 
puted with the median value of the probability distri- 
bution of each parameter. Shaded regions indicate the 
range of models compatible with the 68% confldence in- 
terval for each parameter around the median. In Figure 
8 we show the posteriors of the six torus parameters (the 
vertical shift and the foreground extinction have been 
marginalized) for the galaxy NGC 1097. Those for the 
rest of the Syl galaxies are presented in Appendix A 
(Figures Al to A6). 

The more information the IR SEDs provide, the better 
the probability distributions are constrained. From the 
analysis performed here and in RA09 it appears impor- 
tant to sample the SED in the wavelength range around 
3-4 fim and also at ~18 /Ltm to constrain the model pa- 
rameters. A detailed study of the influence of different 
fllters/wavelengths in the restriction of the clumpy mod- 
els parameter space will be the subject of a forthcoming 
paper (Asensio Ramos et al., in prep.). The posterior 
information for the seven Syl galaxies is summarized in 
Table 8. 

From the individual flts of the Syl galaxies in our sam- 
ple with clumpy torus models we obtain the following 
results: 

1. The average number of clouds along an equatorial 
ray is within the interval A^o=[lj 8], 

^'^ Note that is the foreground extinction from the host 
galaxy, which is different from the Ay^''^ value reported in Ta- 
ble 8, corresponding to the extinction produced by the torus. 



A^CS = 1.086 A^o TV 



mag. 



10 



Ramos Almeida ct al. 



TABLE 7 

Clumpy Model Parameters and Considered Intervals 



Parameter 



Abbreviation Interval 



Width of the angular distribution of clouds a 

Radial extent of the torus Y 

Number of clouds along the radial equatorial direction A^o 

Power-law index of the radial density profile q 

Inclination angle of the torus i 

Optical depth per single cloud ry 



15°, 70°] 
5, 30] 
1, 15] 
0, 3] 
0°, 90°] 
5, 150] 



NGC 109 7 



2 r 



-1 




Av = 2±1 maq 



<i> ^ in 



-1 



NGC 1566 




Wavelength [/v.m] 



Wavelength [/im] 



NGC 6221 



r- 




10 

Wavelength [/J-m] 




10 

Wavelength [/^nn] 



Fig. 6.— High spatial resolution IR SEDs of the Syl galaxies NGC 1097, NGC 1566, NGC 6221, and NGC 6814. Solid and dashed lines 
correspond to the MAP and median models, respectively. Shaded regions indicate the range of models compatible with the 68% confidence 
interval for each parameter around the median. 



2. Low values of a are preferred: a ~ [25°, 45°], 
and intermediate inclination angles of the torus are 
found: i ~ [35°, 85°]. 

3. The radial extent of the torus {Y=Ro/ Rd) is 
weakly constrained within the interval Y ~ [15,20]. 

4. Values in the range from Ty ~ [40,140] are found 
for the optical depth of each cloud for all the galax- 
ies. 

5. The radial density profile appears constrained 
within the interval g=[G.2,1.9], with the only ex- 
ception of NGC 1097, for which q= 



_0 7-1-0.2 



6. The 10 /im silicate feature appears in shallow emis- 
sion or absent in the fitted models with the excep- 
tion of NGC 6221 MAP model. The weak silicate 
feature arises in the clumpy models because both 



illuminated and dark cloud sides contribute to the 
observed spectrum (see Section 5.4 in RA09 for 
a detailed discussion on the silicate feature mod- 
elling) . 

7. The optical extinction produced by the torus along 
the LOS results in Ay^^ < 690 mag for all the 
galaxies. 

8. The foreground extinction from the host galaxy, 
which obscures the AGN direct emission in our 
modelling, results in Ay<5 mag (see Figures 6 and 
7). The values derived are consistent with those 
published in the literature, e.g. the nuclear op- 
tical extinction of Ay =3 mag measured from the 
optical spectrum of NGC 6221 (Levenson et al. 
2001), Ay -1 mag determined from NACO/VLT 
colour maps (Prieto et al. 2005) for NGC 1097, and 
Ay=4.5-4.9 mag reported by MundeU et al. (1995) 
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Fig. 7.— Same as in Figure 6, but for NGC 7469, NGC 3227, and NGC 4151. 




Fig. 8. — Probability distributions resulting from the fit of NGC 1097. Solid and dashed lines represent the mode and the median of each 
distribution and dotted lines indicate the 68% confidence level around the median. The histograms have been smoothed for presentation 
purposes, occasionally leading to small offsets between the solid line corresponding to the mode and the position of the maximum of the 
histograms. 



for NGC 3227. 

All the above intervals or limits of the parameters 
correspond to median values. We chose the medi- 
ans instead of the MAPs because the former gives 
a less biased information about the result, since 
it takes into account degeneracies, while the MAP 



does not. 

The clumpy models succesfuUy reproduce the observed 
Syl SEDs studied here with compatible results among 
them. This is indicating that the NIR and MIR unre- 
solved fluxes employed here are dominated by a combi- 
nation of reprocessed emission from dust in the torus and 
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TABLE 8 

Parameters Derived from the Clumpy Model Fitting 



Galaxy 


cr 

Median 


(°) 

Mode 


Y 

Median 


Mode 


No 
Median 


Mode 


1 

Median 


Mode 


in 

Median 


Mode 


TV 

Median 


Mode 


Median 


(mag) 
Mode 


NGC 1097 
NGC 1566 
NGC 6221 
NGC 6814 
NGC 7469 


29±10 
46±J2 

42±;| 
23±i^ 
41±^ 


28 
52 
36 
20 
38 


19±| 
21±4 
21±5 
18±7 
22±5 


23 
20 
24 
23 
27 


< 2 
2±1 
8±| 
2±1 
3±1 


1 
2 
7 
2 
3 


0.2±0.2 
0.8±g-| 

1.9±0.4 


2.9 
0.1 
0.4 
1.1 
1.9 


36±i7 

59±|| 
43±i 
84±2 (fix) 


43 
86 
79 
39 
84 


42±e9 
104±| 
llOiti 
113±| 
141±B, 


36 
141 
128 
139 
148 


<15 

<690 
<85 
500±120 


1 
160 

1 

1 
430 


NGC 3227 
NGC 4151 


36±|« 

24±r 


33 
19 


19±6 
16±« 


19 
8 


5±3 
2±f 


2 
2 


1.8±0.6 


0.5 
1.8 


49±|| 
43±i 


66 
57 


118±i« 
110±| 


129 
123 


<240 
<60 


1 
1 



Note. — Medians and modes of tfic Syl probability distributions. Those presenting single tails have been characterized with the mode and 
upper/lower limit at 68% confidence. The models include the intrinsic AGN continuum omission. 



direct AGN emission. 

Our modelling results for NGC 1097 somehow contra- 
dict those presented in Mason et al. (2007). The latter 
authors unsuccessfully tried to reproduce the T-ReCS 
11.66 and 18.3 /.tm aperture fluxes that they measured 
for this galaxy using the clumpy models of Nenkova et 
al. (2002). The difference with our result is probably 
due to i) the fact that we obtained unresolved fluxes us- 
ing PSF subtraction over the same images presented in 
Mason et al. (2007), reducing the potential contamina- 
tion from star formation; ii) they did not consider NIR 
data, but only the MIR fluxes; and finally iii) they faced 
the degeneracy problem of the clumpy models without 
using any sophisticated tool as e.g. BayesClumpy. 

5.2.2. Sy2 and intermediate-type Seyfert results 

As a consequence of the publication of the erratum 
Nenkova et al. (2010), where the authors report a mis- 
take in the torus emission calculations, we repeated all 
the fits presented in RA09 using our updated version of 
BayesClumpy. In order to do a proper comparison with 
the results for the Syl galaxies presented here, we per- 
formed the fits of the Sy2 and intermediate-type Seyferts 
considering exactly the same priors as for the Syl. As in 
RA09, we did not consider the direct AGN contribution 
for either Syl.8/1.9 or Sy2. 

For the new fits we considered MIR nuclear fluxes from 
VISIR when available, in addition to the NIR and MIR 
fluxes reported in RA09. We did not include measure- 
ments in the VISIR PAH filters (8.59, 11.25 and 11.88 
/im) for those galaxies with very intense star formation 
such as NGC 7582. Many of the SEDs in RA09 have 
been also updated with NIR data from recent publica- 
tions. In Table Bl (Appendix B) we report the NIR-to- 
MIR SEDs for aU the sample. 

In general, the results from the fitting with the most 
up-to-date version of the interpolated clumpy models, 
which are reported in Table B2 in Appendix B, are com- 
patible with those presented in RA09 at the 1-sigma 
level. Indeed, if we compare the results for the five 
galaxies for which we fitted exactly the same SEDs as in 
the previous work (Circinus, Mrk'573, NGC 1386, NGC 
1808, and NGC 1365) we find that they are practically 
identical. The only fits that are completely different from 
those presented in the previous work correspond to Cen- 
taurus A and NGC 3281. This is a consequence of adding 
new MIR data from VISIR and/or apriori information for 
the inclination angle of the torus (see Appendix B). 



In Table 9 we show the ranges of the parameters found 
for the Syl, intermediate-type Seyferts, and Sy2. We 
have excluded the unrehable fits of NGC 1808 and NGC 
7582 as we did in RA09. In the case of NGC 1808, due 
to the intense star formation that is taking place in its 
nuclear region and to the lower spatial resolution IR SED 
(all from 3-4 m telescopes), it is likely contaminated with 
starlight, resulting in its peculiar shape. For NGC 7582, 
the intense circumnuclear star formation and the edge-on 
orientation of the galaxy make it difficult to isolate the 
torus emission from that of the host galaxy. In the fit, 
the silicate feature is predicted in emission, while from 
MIR spectroscopy shows it in strong absorption (Siebcn- 
morgen et al. 2004). 

We considered the Syl. 8 and Syl. 9 types as a separate 
group in between the Syl & Syl. 5 and Sy2 because the 
IR slopes measured from their SEDs are intermediate be- 
tween those of Syl and Sy2, as measured for our sample 
and also as reported in the literature (Section 4.1 and 
references therein). By looking at the ranges of param- 
eters for the Syl. 8 and Syl. 9 types, we find similarities 
with the Syl & Syl. 5 group in Y, Nq, and Ty, whereas 
a and q are more similar to those of Sy2. 

6. COMPARISON BETWEEN TYPE-1 AND TYPE-2 
SEYFERT NUCLEI. 

The main aim of this work is to enlarge the number 
of Type-1 Seyferts in the original sample of RA09 in 
order to better compare between Type-1 and Type-2 
tori under the assumption that the SEDs studied here 
are torus/AGN dominated. Despite the relatively low 
number of objects considered (7 Type-1 and 9 Type-2 
Seyferts^^), we find that some of the parameters are sig- 
nificantly different between Syl and Sy2. 

To take full advantage of the Bayesian approach em- 
ployed here for the individual fits, the best way to com- 
pare the results for Syl and Sy2 galaxies is to derive joint 
posterior distributions for the full Type-1 and Type-2 
datasets respectively. If Di contains the observed data 
from the i-th SED, assuming that the different SEDs are 
statistically independent, we can use the Bayes theorem 
to write the posterior for all galaxies together as: 

N 

p{e\{D,})^pi{D,}\e)pie) = l[p{D,\9)p{e), (i) 

where 6 ~ (cr, Y, Nq, q, Tv,i)- 

Here we consider the nine Sy2 in RA09 with reliable fits 
(NGC 1808 and NGC 7582 are excluded). 
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TABLE 9 

Ranges of Parameters for SyI, intermediate-type Seyferts, and Sy2 



Type Galaxies cr (°) Y Nq q j (°) Ty 

Syl & Syl.5 7 [25, 45] [15, 20] [1, 8] [0.2. 1.9]^ [35. 85] [40, 140] 

Syl.8 & Syl. 9 3 [35, 70] [20, 25] [1, 7] [0.9, 3.0] [25, 85] [40, 85] 

Sy2 9 [20, 65] [10, 20] [6, 14] [0.0, 3.0] [45, 85] [5, 95] 



Note. — General ranges of parameters resulting from the fits with the clumpy models 
for the seven Type-1 Seyferts, the three Syl.8 and Syl. 9 nuclei, and for the nine Sy2 with 
reliable fits (the extremes of the intervals have been rounded to simplify the comparison). 
All the above intervals correspond to median values. 

''With the exception of NGC 1097, for which g=2.7±{] |. 



Thus, we normahzed all the Syl SEDs at 8.74 /xm and 
fitted them together using BayesClumpy, and we did the 
same for the Sy2. For those galaxies without flux mea- 
surements in the Si-2 filter we performed a quadratic 
interpolation of the SED and used the interpolated val- 
ues at 8.74 /.tm to normalize the real data. We did not 
use the interpolated values in the fits. We considered the 
mean redshift for the Syl (z=0.0061±0.0045) and for the 
Sy2 (0.0078±0.0051) in the fits^^. In Figure 9 we show 
the Syl (left panel) and Sy2 fits (right panel). Note that 
the IVIAF and median models predict a flat SED with the 
silicate feature in weak emission for the Syl galaxies, and 
steeper and with the silicate band in shallow absorption 
for Sy2. 

The comparison between the Syl and Sy2 posterior 
distributions is shown in Figure 10. From a visual inspec- 
tion it is clear that the joint posteriors of the parameters 
No, q, TV, and a are completely different between Syl 
and Sy2. There is not overlap between the 1-sigma in- 
tervals. In Table 10 we report the median and mode 
values of the histograms in Figure 10. 

In order to quantify how different the probability dis- 
tributions are, we calculated the KuUback-Leibler diver- 
gence (KLD; Kullback & Leibler 1951) between the Syl 
and Sy2 posteriors. This divergence takes into account 
the full shape of the posterior and it is always a pos- 
itive value, and it is equal to zero when two distribu- 
tions are identical. Therefore, the larger the value of 
KLD, the more different the posteriors. We find KLD>1 
for a (KLD=5.2), A^o (KLD=25), q (KLD=1.2), and Ty 
(KLD=21). These are indeed the four parameters whose 
1-sigma regions do not overlap (see Figure 10) and thus 
we consider their differences significant between Syl and 
Sy2^^. For both Y and i we find KLD<1 and similar 
median values between Syl and Sy2. 

Syl tori are narrower and have fewer clouds 

(cr=44°±8°; No=4±l) than those of Sy2 {(j^m°±^°o] 

No=ll±^). The radial density distribution of the clouds 
is also different between the two Seyfert types accord- 
ing to this analysis: in Sy2, the majority of the clumps 
are distributed very close to the nucleus (i.e. steep ra- 
dial density distribution; g=2.3±0.1) whereas for Syl the 
clouds distribution is flatter (g=:0.8±0.2). On the other 
hand, the optical depth of the clouds in Syl tori is larger 
(ry = 133±8) than in Sy2 (Ty=30±l). 
By taking a closer look to the right panel of Figure 9, 

Indeed, since all the galaxies are local Seyferts, the results arc 
the same if we consider that the SEDs are rest-frame 

^"^ See further discussion on the Sy2 q parameter results below. 



which corresponds to the Sy2 fit, it is clear that some 
of the data points are underestimated by the model. 
This happens because of the Circinus SED, which has 
the shortest errors bars, has more weight in the fit than 
the rest of the SEDs. In order to check that our Sy2 
results are not completely biased by Circinus, we have 
repeated the fit excluding it, as well as NGC 1386, which 
is the least restricted SED in terms of data points, get- 
ting rid of the extremes. From the new fit, we find even 
larger differences with the Syl values of cr, A^o, and ry. 
On the other hand, the q joint posterior becomes com- 
parable to that of the Syl, probably as a consequence 
of getting rid of the two SEDs fitted with the largest q 
values among the Sy2. Considering all the previous we 
prefer to be cautious about the q parameter, and only 
consider cr, iVg, and ry genuinely different between Syl 
and Sy2. 

Interestingly, we find high as well as low values of the 
inclination angle of the torus for Syl and Sy2 (see Table 
9). This variety in the i values translates into the similar 
median values found for the joint Syl and Sy2 poste- 
rior distributions (47°±^o for Syl and M°±f°o for Sy2), 
which are also intermediate within the considered prior 
(i=[0°,90°]). This is telling us that, in the clumpy torus 
scenario, the classification of a Seyfert galaxy as a Type- 
1 or Type-2 depends more on the intrinsic properties of 
the torus rather than in its inclination. 

Our results contradict those presented by Honig ct al. 
(2010), who find similar averaged values of Nq for both 
the Syl and Sy2 galaxies in their sample. However, it 
is worth mentioning that they fixed low values of the 
inclination angle of the torus for Syl and large for Sy2, 
what can likely have infiuenced their results. 

In Figure 11 we represent the median values of cr 
and No for the different Seyfert types over the cover- 
ing factor contours-'^^. The covering factor is defined as 
Ct = 1 - / e-^^o-s(*)dcos(i). Type-1 nuclei tend to be 
located within lower Ct contours (Ct <0.6) than those 
of Type-2s, for which Ct >0.5, with the exception of 
Centaurus A and Mrk 573^^. We have represented with 

A similar plot showing values for Sy2 galaxies from RA09 
and PG quasars from Mor et al. (2009) was shown in the talk by 
M. Elitzur at the Physics of Galactic Nuclei conference held 15-19 
June, 2009 at Ringberg Castle. Proceedings published online at 
http: / /www. mpe.mpg.de/events/pgn09/online_ proceedings.html. 

As discussed in RA09, the fit of Centaurus A is complicated 
by the presence of a dust lane of Ay ~7-8 mag that is likely af- 
fecting the NIR nuclear fluxes, as well as the possible synchrotron 
contamination of the MIR fiuxes. Mrk 573 has been recently re- 
classified as an obscured Narrow-line Seyfert 1 (NLSyl) based on 
NIR, spectroscopy (Ramos Almeida et al. 2008, 2009b). However, 
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Fig. 9. — Same as in Figure 6, but for the Syl (left) and Sy2 SEDs (right) normalized at 8.74 /im. 



TABLE 10 

Statistics of the comparison between SyI and Sy2 parameters 



Type 


cr 
Syl 


(°) 
Sy2 


Y 

Syl Sy2 


Syl 


Sy2 


Syl Sy2 


i 

Syl 


(°) 

Sy2 


TV 

Syl 


Sy2 


Medians 
Modes 


44±8 
42 


63±| 
65 


21±4 23±| 
21 29 


4±1 
4 


11±2 
11 


0.8±0.2 2.3±0.1 
0.9 2.2 


47±| 
46 


54±}; 
55 


133±| 
132 


30±1 
30 



larger symbols the median values from the joint a and 
iVo posteriors reported in Table 10. 

Since the covering factor is a non-linear function of the 
torus parameters, we took full advantage of our Bayesian 
approach and generated joint posterior distributions for 
Ct from those in Figure 10, which are shown in the 
left panel of Figure 12. The median values of the his- 
tograms are CT(Sy2)==0.95±0.02 and CT(Syl)=0.5±0.1. 
The divergence between the Syl and Sy2 Ct posteriors 
is KLD=28, indicating that the difference is significant 
(the 1-sigma regions do not overlap). Thus, Syl tori in 
our sample have lower CtS than those of Sy2, implying 
that they arc intrinsically different. 

6.1. Near- infrared Emission and Torus Angular Width 

As reported in Section 4.1, Type-1 Seyferts present 
characteristically higher H/N ratios and flatter SEDs 
than those of Sy2 nuclei. The enhancement on the NIR 
emission of Type-1 AGN is produced by the hot dust 
from the directly-illuminated faces of the clumps in the 
torus which are close to the central engine, and also by di- 
rect AGN emission (i.e., the tail of the optical/ultraviolet 
power-law continuum), which strongly flattens their IR 
SEDs (Rieke & Lebofsky 1981; Barvainis 1987). 

Based on the IR SEDs presented here, the relative NIR 
contribution to the SED is generally correlated with the 
Seyfert type. Sy2 galaxies show lower NIR to MIR ratios 
(H/N=0.003±0.002) than Syl (0.06±0.03). 

In the context of the clumpy models, the presence of a 
cloud along the LOS, which may occur from any viewing 
angle, results in a Type-2 classification. Cloud encoun- 
ters are more probable at large inclination angles (i), but 
there is always a finite probability for having an unob- 
scured view of the AGN. In fact, the likelihood of inter- 
here we considered it as a Sy2 because of the similarity in SED 
shape with the rest of the sample. 



cepting a cloud along a LOS depends on the combina- 
tion of i, Nq, and a. Thus, the preference for lower val- 
ues of these parameters (especially A^o and a) in Type-1 
Seyferts increases the likelihood of unimpeded views of 
some directly- illuminated cloud faces (i.e., those on the 
"back" side of the torus) and direct AGN emission, re- 
sulting in an increase of the NIR flux. The latter can be 
represented in terms of the escape probability Pesc (see 
equation 4 in Nenkova et al. 2008a). For a total number 
of clouds Ntos along a path, Pesc — exp{~NLos) when 
the clouds are optically thick {t\ > 1). 

In Figure 13 we show the dependence of the H/N ratio 
on the escape probability. All the Sy2 are in the bottom- 
left corner of the diagram, whereas the Syl have higher 
values of the H/N ratio and Pesc=[l%,92%]. Syl.8 and 
Syl. 9 galaxies present intermediate values between those 
of Syl and Sy2. The derived joint posterior distributions 
of the escape probabilities for Syl and Sy2 (KLD=29 
between them) arc shown in Figure 12, and the median 
values of the histograms are Pesc(Syl) = 18±3% and 
Pesc(Sy2) ^ 0.05±0:P%. 

Thus, while for tori with high values of i, A^o- and a the 
probability of having a direct view of the AGN is very 
little, that increases for objects with narrower and less 
inclined tori, and containing less clumps. In this work we 
show for the first time that, in the clumpy torus scenario, 
the classification as a Type-1 or Type-2 Seyfert depends 
more on the intrinsic properties of the torus, rather than 
in the inclination angle itself. The Syl galaxies in our 
sample have larger Pesc than the Type-2 Seyferts, and as 
a consequence of that, we detect the broad lines in their 
spectra. 

6.2. AGN Luminosities 

The clumpy model fits yield the intrinsic bolometric lu- 
minosity of AGN (L^^^) by means of the vertical shift 



Testing the Unification Model in the Infrared 



15 




Fig. 10. — Same as in Figure 8, but for the joint Syl and Sy2 SEDs. KLD values derived from the comparison between Syl and Sy2 for 
each parameter are labelled. 



applied to match the observational data points. Combin- 
ing this value with the torus luminosity (L^°[), obtained 
by integrating the corresponding model torus emission 
(without the AGN contribution), we derive the repro- 
cessing efficiency (RE) of the torus (L^°[/Lf„"f^). The 
previous values are calculated on the Bayesian frame- 
work, by combining the posterior distributions of the 
model parameters. Median values and 1-sigma intervals 
for Syl galaxies are reported in Table 11, and those for 
Sy2 and intermediate-type Seyferts are shown in Table 
B3 (Appendix B). 

Sy2 tori in our sample are more efficient reprocessors 
than Syl, absorbing and re-emitting the majority of the 
intrinsic AGN luminosity in the IR: RE(Sy2)=[0.4, 1.0], 
with a median value of 0.8 and RE(Syl)=[0.2,0.7], with 



median of 0.5. It makes no sense to compare the derived 
joint Syl and Sy2 posterior distributions of RE, because 
we normalized the SEDs to perform the fits, and conse- 
quently the derived L^^^ arc meaningless. 

We considered a possible dependency of the RE (or al- 
ternatively the covering factor; Ct) on L^'^^, since the 
amount of incoming radiation from the AGN could pos- 
sibly have some infiuence on the reprocessed energy or 
even in the torus properties (e.g., receeding torus sce- 
nario; Lawrence et al. 1991). However, we find no re- 
lationship between the two quantities in the luminosity 
range considered. This means that the reprocessing effi- 
ciency depends primarily on the total number of clouds 
available to absorb the incident radiation, i.e. on the 
torus covering factor. However, the possible dependence 
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Fig. 11. — a versus No for the individual galaxies. Either median values or upper/lower hmits are taken from the fits presented here. 
Dots eorrespond to Type-1 Seyferts, triangles to Syl.8 and Syl.9, and squares to Sy2. Error bars indieate 68% eonfidence level around 
the median. Note the segregation between Seyfert types, indieating the intrinsic difference between their tori in terms of covering factor 
(indicated in contours). The big dot and square correspond to the average a and Nq values for Syl and Sy2 from Table 10. 
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Fig. 12. — Joint posterior distributions of the torus covering factor (left panels) and escape probability (right panels) for Syl (top) and 
Sy2 galaxies (bottom). The values of the KuUbaek-leibler diyergcncc obtained from the comparison between Syl and Sy2 are KLD=28 for 

(3i|^tfisitHIiB=42igQ|^eifliQ§. on the AGN luminosity consider- ing a broader luminosity range is further investigated in 
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Fig. 13. — H/N versus Peso- Lower values of A^o, cr, and i result in higher Pesc, resulting in higher NIR-to-MIR ratios. The inset shows 
an amplification of the region occupied by the Sy2. Symbols are the same as in Figure 11. The big dot and square correspond to the 
median values of Pesc for Syl and Sy2 from Figure 12, and to the H/N values of the Syl and Sy2 average templates (Table 6). 



TABLE 11 
BoLOMETRic Luminosity Predictions 



Galaxy 



bol 



hoL I bol 



Ro (pc) 



T AGN 



NGC 1097 
NGC 1566 
NGC 6221 
NGC 6814 
NGC 7469 



2 4_|_0.8 

6-2±|:o 
8.1±|« 
3.7±0.8 



NGC 3227 
NGC 4151 



1.3±, 



^■"^0.9 



X 


10*2 


0.4±0.1 




1.0 


X 


10^2 


0.4 


a 


X 


10*2 


0.4±0.1 


0.6±0.2 


6.3 


X 


10*2 


1.1 


b 


X 


1042 


n 7+0.5 
'^■'='=0.3 


n 7-I-0-3 


1.3 


X 


10« 


2.1 


c 


X 


10*2 


0.4±0.2 


0.7±0.3 


2.7 


X 


10« 


3.3 


d 


X 


Iq44 


0.5±0.1 


5±i 


3.4 


X 




0.9 


e 


X 


10« 


n 7-1-0. b 


0.9±h 


3.8 


X 


lo-i^ 


2.9 


f 


X 


10« 


0.5±0.1 


1 4-I-0-8 


1.7 


X 


Iq44 


3.7 


g 



References. — (a) Terashima et al. (2002); (b) Levenson et al. (2009); (c) Levenson et 
al. (2001); (d) Gandhi et al. (2009); (e) Nandra et al. (2007); (f) Lamer et al. (2003); (g) 
Bcckmann et al. (2005). 

Note. — Bolometric luminosities corresponding to the AGN luminosity (Lj^ 



AGN. 



erg s~^). Columns 3 and 4 correspond to the RE (L^°[/ijJ^^'^) and the outer radius of 
the torus calculated using LjJ^^'^ and Y . Absorption-corrected 2-10 keV X-ray luminosities 
are taken from the literature (references below). L^^j^ (erg s~^) is derived from 20xL^'''^. 



Alonso-Herrero et al. (2011). 

Figure 14 shows that there is a correlation between RE 
and the torus covering factor for the galaxies in our sam- 
ple. The larger Ct the more efficient reprocessor. Type-2 
tori have in general larger RE than those of Type-1, with 
the exceptions of Centaurus A and Mrk 573. Despite the 
relatively small sample, Figure 14 shows a segregation 
between Syl and Sy2 galaxies in terms of reprocessing 
efficiency and covering factor. 

// all Seyfert nuclei are identical, as the unified model 
predicts, only the viewing angle should determine the 
classification, not the properties of the torus itself. While 
our results are limited by the small sample analyzed, they 
suggest instead that the Type-1 /Type-2 classification de- 
pends on the torus intrinsic properties rather than in the 
mere torus inclination. 



The bolometric luminosity of the intrinsic AGN de- 
rived from the fits (L^^^; column 2 in Table 11) can 
be directly compared with those from the absorption- 
corrected 2-10 keV luminosities compiled from the litera- 
ture {Lfi;^ ; column 6 in Table 11), which is an effective 
proxy for the AGN bolometric luminosity (Elvis et al. 
1994). To obtain L^f J from the intrinsic 2-10 keV val- 
ues we applied a bolometric correction factor of 20 (Elvis 
et al. 1994). 

We find similar values of L^^j)^ and L^^^ for all the 
Syl (sec Table 11). In the case of NGC 1097, we expect to 
have some contamination from the nuclear starburst (see 
Section 5.2.1 and Mason et al. 2007). Thus, if the L^^j)" 
value represents the AGN contribution-only, Lf^^ may 
be overestimated because of the starburst contribution, 
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Fig. 14. — Reprocessing efficiency versus torus covering factor for the individual galaxies. Higlier values of Ct translate into more efficient 
reprocessors. In general, Type-2 tori are more efficient tfian Type-1 tori, with the exceptions of Ccntaurus A and Mrk 573. Symbols are 
the same as in Figure 11. 

thus resulting in a low ratio between the two measure- 
ments. The good agreement with the observational X-ray 
measurements reinforces the results from our SED mod- 
elling. The same comparison with X-ray data for Sy2 
and intermediate-type Seyferts is shown in Table B3 in 
Appendix B. 

6.3. Torus Size 

In general, the IR SED fitting does not constrain the 
size of the torus (Y) as well as other model parameters 
(see Section 5.3 in RA09). The NIR and MIR obser- 
vations are sensitive to the warm dust (located within 
^10 pc of the nucleus), which depends on the combina- 
tion of model parameters Nq, q, and Y (Thompson ct al. 
2009). FIR observations would be more sensitive to the 
torus extent independently. However, the fits performed 
here for Type-1 and Type-2 Seyferts arc consistent with 
a small torus size, confined to scales of less than 6 pc (see 
below) . 

Uniform density models require the dusty torus to ex- 
tend over large dimensions, to provide cool dust that pro- 
duces the IR emission (e.g., Granato & Danese 1994). In 
contrast, in a clumpy distribution, different dust temper- 
atures can coexist at the same distance, including cool 
dust at small radii (Nenkova et al. 2002), so large tori, 
which are inconsistent with imaging and interferometry 
results, are not necessary. Indeed, for the Seyfert galax- 
ies considered here, we found Y ranging from 10 to 25 
(see Table 8), showing that small tori can account for the 
observed IR nuclear emission. 

The outer size of the torus scales with the AGN 
bolometric luminosity: Ro = YRd, so assuming 
a dust sublimation temperature of 1500 K, Ro = 
0.4 Y (L^^^ /IQ-^^f-^ pc. We derived Ro posterior dis- 
tributions from those of L^^^ and Y and find that all 
tori in our sample have outer radii smaller than 6 pc (Ta- 



bles 11 and B3), in agreement with MIR direct imaging 
of nearby Seyferts (Packham et al. 2005; Radomski et 
al. 2008) and also interfcrometric observations (Jaffe ct 
al. 2004; Tristram et al. 2007; Meisenheimer et al. 2007; 
Raban et al. 2009). 

For example, the estimated outer radius for NGC 1097 
(Ro=0.4±Q-2 pc) is in agreement with the value derived 
from NIR NACO/VLT observations (Prieto et al. 2005), 
which placed the radius of the central compact source in 
r<5 pc. Mason et al. (2007) also derived an upper limir of 
19 pc for the size of the unresolved component from the 
MIR images employed in this work for NGC 1097. For 
the case of NGC 7469, for which we derive Ro=5±2 pc, 
Tristram et al. (2009) reported an estimation of the size 
of the dust distribution of 10 pc, based on MIDI/VLT 
interfcrometric observations, although compromised by 
the high level of noise and the lack of fringes. 

7. CONCLUSIONS 

We present new subarcsecond resolution MIR imaging 
data at 8.7 and 18.3 /im for the three Type-1 Seyfert 
galaxies NGC 6221, NGC 6814, and NGC 7469. NGC 
7469 and NGC 6221 appear extended, with part of this 
extended emission associated with emitting-dust heated 
by star formation. On the contrary, NGC 6814 lacks of 
any extended emission. Nuclear MIR and NIR fluxes for 
the three galaxies as well as for NGC 1566, NGC 1097, 
NGC 3227, and NGC 4151 arc reported. We construct 
nuclear SEDs that the AGN dominates and fit them with 
clumpy torus models and a Bayesian approach to derive 
torus parameters. The main results of this work for the 
individual Syl galaxies and from the comparison with the 
Sy2 and intermediate-type Seyferts in RA09 arc summa- 
rized as follows: 

• We derived an average Type-1 Seyfert template 
from the individual Syl SEDs, which is flatter 
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(mean IR slope am = 1.7 ± 0.3) than the Type-2 
mean SED presented in RA09 {am = 3.1 ± 0.9). 

• The NIR-to-MIR flux ratios measured from the in- 
dividual SEDs are larger for Syl (0.07±0.03) than 
for Sy2 (0.003±0.002). Indeed, the distributions of 
NIR-to-MIR values are significantly different be- 
tween the two types at the 100% confidence level. 

• The interpolated version of the clumpy models of 
Nenkova et al. (2008a, b) employed here success- 
fully reproduces the nuclear IR SEDs of Type- 
1 Seyfcrts with compatible results among them. 
Consequently, the observed nuclear IR emission of 
these galaxies can be accounted for by dust heated 
by the central engine and direct AGN emission. 

• We derive joint posterior distributions for Syl and 
Sy2 and find that the differences in Ng, ry, and a 
between Type-1 and Type-2 tori are significant ac- 
cording to the Kullback-Leibler divergence and the 
lack of overlap between their 1-sigma confidence 
intervals. 

• We find that Syl tori are narrower and have fewer 
clouds than those of Sy2. Additionally, the optical 
depth of the clouds in Syl tori is larger than in Sy2. 

• There is not a clear trend in the values of the incli- 
nation angle of the torus for Syl and Sy2 (slightly 
larger values are found for Sy2). 

• The larger the covering factor of the torus, the 
smaller the likelihood of intercepting a cloud along 



a LOS. In our sample, Seyfert 2 tori have larger cov- 
ering factors and smaller escape probabilities than 
those of Seyfert 1. 

• Despite the limited number of galaxies considered, 
we find that Type-2 tori are in general more effi- 
cient reprocessors than those of Type-1. Indeed, 
there is a correlation between the reprocessing effi- 
ciency and the torus covering factor. 

• For the Seyfert galaxies studied here, we find that 
tori with outer radii smaller than 6 pc can account 
for the observed NIR/MIR nuclear emission, in 
agreement with MIR interferometric observations. 

Summarizing, we find tantalizing evidence, albeit for a 
small sample of Seyfert galaxies and under the clumpy 
torus hypotesis, that the classification as a Type-1 and 
Type-2 depends more on the intrinsic properties of the 
torus than on its mere inclination towards us. 



APPENDIX 

FITTING RESULTS FOR SYl GALAXIES 

Here we include the posterior distributions of the Syl galaxies NGC 1566, NGC 6221, NGC 6814, NGC 7469, NGC 
3227, and NGC 4151 (Figures Al to A6). 
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TABLE Al 

VISIR MIR FLUXES FROM THE LITERATURE FOR THE SyI GALAXIES 



Galaxy 


Flux Density (mjy) 


Filtor(s) 


Reference(s) 


NGC 1097 


28±7, 49±12 


12.27, 18.72 


B1,B2 


NOG 1566 


63±9, 128±32 


11.88, 18.72 


B2 


NGG 6814 


99±6, 96±6 


11.25,13.04 


B3 


NGG 7469 


448±16. 595±18, 630±17, 1270±317 


10.49,12.27,13.04,18.72 


B1,B4,B2 



NGG 3227 180±11, 320±22 8.99, 11.88 B4 



References. — (Bl) Horst ct al. (2008); (B2) Rcunanen et al. (2010); {B3) Gandhi et 
al. (2009); (B4) Honig ct al. (2010) 
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Fig. A4. — Same as in Figure 8, but for the galaxy NGC 7469. 




Fig. A5. — Same as in Figure 8, but for the galaxy NGC 3227. 
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TABLE Bl 

NIR AND MIR FLUXES FROM THE LITERATURE FOR Sy2 AND INTERMEDIATE-TYPE SeYFERTS 



Galaxy 


Seyfert 


Flux Density (mjy) 


Filter(s) 


R,cfcrence(s) 


CenA 


2 


1 Q-l-n 1 d ^-l-n '44-1-9 900-1-40 


NAPO T H K T,' 


Al 






^4*^-1-97 Q47-t-9Q 1100-1-100 14^1-1-7*^ 9*^00-1-^17^1 


\/"Tc;t'r 1 n 4q 1 1 9^^ 1 1 rj^ 1 9 97 1 s 79 


Bl B2 


r<- ■ 


2 


1 fi-t-0 9 4 S-l-0 7 1Q-I-9 '^l-t-'^ '^SO-I-'^S 1Q00-I-1Q0 




A2 


TP f^Ofi'^ 
ouuo 


2 


n "^-1-0 1 4 77-1-0 


TTI fiflW 'Fr999M 


A "4 A 1 






finQ-l-99 797-1-9^^ Q9^-l-9^ lO'^fi-i-^fi 


VmTR 1 4Q 1 T 9"^ 1 1 1 9 97 


133 


iVl 1 iV O 1 ij 


2 


n 1 c;-|-0 OR ^4-1-0 04 9-1-0 R IS S-t-"^ S 41 "^-l-S 


FIIOWFIfiOW NSFCam K' T, M 


A4 


NGC 1386 


2 


0.2±0.1 


F160W 


A3 




2 


15 5-1-4 5 30 5-1-8 5 27 5-1-10 5 


ISAAC J,Ks,L' 


A5 




2 


22-1-0 1 3 


pi ROW 


A3 






138±11 


VTSTR 1 3 04 


B4 


NGC S2S1 


2 


1 '5-1-0 2 7 7-1-0 8 103-1-9 207-1-25 


TRAG-1 HK TRCAM3 T,' M 


A6 






481±24, 1016±52 


VISIR 11.25,13.04 


B4 


NGC 4388 


2 


0.06±0.02, 0.71±0.28, 40±8 


F110W,F160W, NSFCam L 


A4 






147±11, 375±28 


VISIR 11.25,13.04 


B4 


NGC 7172 


2 


<0.4, 3.4±0.7, 30±6, 61±12 


IRCAM3 H,K,L',M 


A7 






165±27 


VISIR 12.27 


Bl 


NGC 7582 


2 


11±1, 18±2, 96±10, llOitll, 142±21 


F160W, NACO 2.06,L,4.05, ISAAC M 


B2,A8 






236±27, 550±130 


VISIR 8.99, 18.72 


B3,B2 


NGC 1365 


1.8 


8.3±0.8 


F160W 


A9 


NGC 2992 


1.9 


<1, 2.8±0.6. 22.7±4.5, 35.7±7.1 


IRCAM3 H,K,L',M 


A7 






312±31 


VISIR 11.25 


B5 


NGC 5506 


1.9 


13±1, 53±1, 80±1, 290±10, 351±29 


NACO J,H,K,L', NSFCam M 


B2.A11 






900±100, 1400±350 


VISIR 11.88, 18.72 


B2 



References. — (Al) Mcisonhoimor ot aL (2007); (A2) Prioto ot al. (2004); (A3) Quillon ot al. (2001); (A4) Alonso-Herroro ot al. (2003); 
(A5) Galliano & AUoin (2008); (A6) Simpson (1998); (A7) Alonso-Herrero et al. (2001); (A8) Prioto ot al. (2002); (A9) Carollo ot al. (2002); 
(AlO) Kulkarni et al. (1998); (All) Ward ot al. (1987); (Bl) Horst ot al. (2008); (B2) Prioto ot al. (2010); (B3) Honig ot al. (2010); (B4) 
Gandhi et al. (2009); (B5) Haas ot al. (2007) 

Note. — Ground-basod instruments and telescopes are: NACO and ISAAC on the 8 m VLT, NSFCam on the 3 m NASA IRTF, IRCAM3 
on the 3.8 m UKIRT, and IRAC-1 on the 2.2 m ESQ telescope. Measurements in the FllOW, F160W, and F222M filters arc from NICMOS 
on HST. 



TABLE B2 

Parameters from the fits for Sy2 and intermediate-type Seyferts 



Galaxy 


a Y No 


9 


i 


TV 






Med Mode Med Mode Med Mode 


Med Mode 


Med Mode 


Med Mode 


Med Mode 



CenA 
Circinus 
IC 5063 
Mrk 573 
NGC 1386 
NGC 1808 
NGC 3081 
NGC 3281 
NGC 4388 
NGC 7172 
NGC 7582 



20±2 

63±4 
51±10 
30±18 
49±}i 
26±P 
55±°7 
34±3 
60±'5 
58± 
34± 



35±f" 

>68 



20 

68 
56 
23 
52 
20 
61 
34 
67 
66 
23 



13±2 

12±| 
17±8 
17±8 
18±^ 
18±| 
10±2 
20±5 

14±| 
18±|_ 



13 
27 
10 
25 
17 
18 
22 
9 
22 
8 
18 



13±^ 
11±2 

13±2 

11±2 

8±4 

10±3 

14±i 

11±2 

10±f 
3±1 



15 
11 
14 

5 

13 
12 
12 
15 
14 
9 
2 



0.4±0 

2.6±0 

0.8±« 

0.9±« 

1.5± 

1.3±« 

1.4±0 



<0.2 
0.5±0'5 
1.9±' 



2.3± 



0.5 
2.7 
0.3 
0.6 
2.3 
0.9 
2.3 
0.0 
0.2 
2.5 
2.6 



70±2 
85±2 (fix) 
81±| 

85±2 (fix) 
85±2 (fix) 
41±P 
54±2i 

63±| (a) 

fiQ-LlB 

43± 



69 
85 
85 
85 
85 
53 
79 
60 
28 
77 
45 



45±if 
31±2 
90±19 
72±« 

lll±|f 

29±?| 
15±2 
31±16 
7±? 



i6±;. 



86±|S 
39±| 
42±3 



42 

31 
96 
53 
51 
122 
21 
14 
25 
6 

13 



235diQQ 
365±™ 

1 1 or:_|_315 

805±44Q 

<100 

<230 

120±20 

<270 

50±15 

<15 



200 
340 
990 
105 
425 
1 

70 
120 
110 
50 
1 



NGC 1365 
NGC 2992 
NGC 5506 



32 
28 



18±7 22 
18±7 24 
23±| 26 



7±4 
<2 



0.9±i 
>2.9 



0.7 
0.2 
3.0 



27±; 

85±2 (fix) 



19 
78 
85 



96 
34 
43 



<110 

2io±;« 

60±8 



1 

235 
60 



Note. — For the galaxies Circinus, Mrk 573, NGC 1386, and NGC 5506 tlie i parameter has been introduced as a Gaussian prior into the 
computations, centered at 85° with a width of 2°, based on other observations. Probability distributions presenting a single tail have been 
characterized with the mode and upper/lower limits at 68% confidence. Syl.8 and Syl.9 have been fitted with the geometry corresponding to 
torus emission-only. 

''The inclination angle of the torus was introduced as a uniform prior j=[60°,90°] following Storchi-Bcrgmann et al. (1992). 



FITTING RESULTS FOR SY2 AND INTERMEDIATE-TYPE SEYFERT GALAXIES 

Here we include the results of the fits of the intermediate- type Seyferts and the Sy2 galaxies in RA09 (Figures Bl 
and B2) using the data reported in Table Bl and the most up-to-date version of the Nenkova et al. (2008a, b) models 
(see erratum Nenkova et al. 2010). The resulting model parameters for the individual galaxies are shown in Tables B2 
and B3. In Figure B3 we show the posterior distributions from the fit of Centaurus A (see the electronic edition of 
the Journal for the rest of the Sy2 and intermediate- type Seyferts posteriors.). 

As mentioned in Section 5.2.2, the results presented in this Appendix are compatible in general with those presented 
in RA09 at the 1-sigma level. For the majority of the galaxies, we have included IR data from recent publications. 
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Fig. B1.— Same as in Figure 6, but for the Sy2 galaxies Centaurus A, Circinus, IC 5063, Mrk 573, NGC 1386, NGC 1808, NGC 3081, 
and NGC 3281. 
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1 10 1 10 

Wavelength [/J-m] Wavelength [jJ-ru] 

Fig. B2.— Same as in Figure 6, but for the Sy2 galaxies NGC 4388, NGC 7172, and NGC 7582; the Syl.8 NGC 1365, and the Syl.9 
galaxies NGC 2992 and NGC 5506. 

In particular, the use of VISIR MIR data have been key in predicting the silicate feature in absorption in the case of 
Ccntaurus A^o and NGC 3281, and in emission in NGC 3227. 

In the case of Ccntaurus A. the fit presented in RA09 predicted the silicate feature in emission, in clear contradiction 
with MIR spectroscopic data (Siebenmorgen et al. 2004; Meisenheimer et al. 2007). By including the VISIR data for 
this galaxy, the fitted models predict a broad silicate feature, in agreement with the observations. Additionally, we 
have also introduced the inclination angle of the torus as a gaussian prior centred in 70° as infered from the inclination 
of the inner radio jet of Ccntaurus A (Tingay ct al. 1998). This value is also coincident with the 66° inclination angle 
of the obscuring structure derived from MIR interferometric observations (Burtscher et al. 2010). 

For the Sy2 galaxy NGC 3281 we also have considered a uniform prior in the range [60°, 90°] for the inclination 
angle of the torus, based on the inclination of the ionization cone axis derived in Storchi-Bcrgmann et al. (1992) from 
optical imaging and spectroscopic observations, which implies a nearly edge-on orientation of the obscuring structure. 



'^^ In order to account for the extinction caused by the dust lane of Ay ~ 7-8 mag in Ccntaurus A, we have considered the Chiar & 
Tielens (2006) law in the fit fixing Ay=S mag. 
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0.5 1.0 1.5 20 30 dO 50 60 70 8S 68 70 72 71 

q T, 1 

Fig. B3. — Same as in Figure 8, but for the galaxy Centaurus A. See the electronic edition of the Journal for the rest of the Sy2 and 
intermediate- type Seyfcrts. 



TABLE B3 

BoLOMETRic Luminosity Predictions 



Galaxy 



bol 



T AGN 
'^Xbol 



Xhol ' bol 



Centaurus A 


4.7± 


Circinus 


l.Oib 


IC 5063 


2.7± 


Mrk 573 


1.5± 


NGC 1386 


2.8± 


NGC 1808 


3.1± 


NGC 3081 


7.6± 


NGC 3281 


1.1± 


NGC 4388 


2.3± 


NGC 7172 


8.4± 


NGC 7582 


1.4± 


NGC 1365 


8.1± 


NGC 2992 


3.0± 


NGC 5506 


9.3± 



X 10^2 

X 10« 
X lO** 
X lO-** 

X 10*2 
X 10-*2 
X 10-*2 

X lO** 
X 10« 

X 10*2 

X 10« 



0.4±0.1 
0.8±0.1 

0.4±0.2 

0.8±»:| 
0.9±0.3 
0.7±0.1 
1.0±0.3 

0.6±0.1 



0.4±0.1 

o.8±0:3 

^•'-'='=0.8 

2.8±J-| 
0.4±0'.2 
0.4±0.2 
0.7±0.3 
1.3±0.3 
1.3±0.4 

0.8±0.4 



X 10*2 
X 10« 
X 10« 



0.6±o| 
0.5±0.1 



1.3 X 10*=^ 


3 


a 


1.2 X 10*3 


1.2 


b 


1.7 X 10** 


0.6 


c 


4.4 X 10** 


3 


d 


1.3 X 10*3 


5 


e 


2.2 X 10*1 


0.1 


f 


1.0 X 10** 


13 


d 


3.0 X 10** 


3 


g 


1.5 X 10** 


6 


h 


1.1 X 10** 


13 


i 


9.8 X 10*3 


7 


j 


3.0 X 10*3 


4 


k 


3.0 X 10*3 


1.0 


1 


2.2 X 10** 


2 


m 



0.7±0.3 

2.9±»:| 



References. — (a) Markowitz ct al. (2007); (b) Soldi ct al. (2005); (c) Turner ot al. (1997); (d) 
RA09; (e) Levonson et al. (2006); (f) Jimcncz-Bailon ct al. (2005); (g) Vignali &; Comastri (2002); 
(h) Elvis ct al. (2004); (i) Awaki ct al. (2006); (j) Turner ct al. (2000); (k) Risaliti et al. (2005); (1) 
Yaqoob ct al. (2007); (m) Lamer ct al. (2000); (n) Lamer et al. (2003); (o) Beekmann et al. (2005). 

Note. — Same as in Table 11, but for the Sy2 and intermediate-type Seyfert galaxies in RA09. 
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